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NOMENCLATURE 


Bore  radius  of  a mandrel,  present  design  a = 0.9845". 
Outside  radius  of  a mandrel,  present  design  b = 1.795". 
Radius  of  a semi-circular  crack;  also  the  depth  of  a 
frontal  crack. 

Modulus  of  elasticity  in  tension  and  in  shear  respectively. 
Uniform  radial  stress  transferred  to  the  contact  region  of 
the  mandrel  h = 38  Ksi  used  in  this  report  to  non-dimen- 
sionalize  other  stresses. 

Thermal  conductivity  of  steel. 

Stress  intensity  factors  of  mode  I. 

Cylindrical  coordinate. 

Ratio  of  b/a. 

Compressive  stress  in  the  x-direction. 

Temperature,  Ta  temperature  at  r = a. 

Displacement. 

(l-2v)/(2-2v). 

Coefficients  in  linear  approximation  of  tangential  stress. 
One-half  of  the  angle  of  contact  between  the  mandrel  and 
the  tube. 

Poisson's  ratio. 

Ratio  of  a/b,  p = 1/R,  for  the  present  design  p ~ 0.55. 
Normal  stress. 


Shear  stress. 


n = 0,  1,  ....n,  stress  functions. 

Partial  derivatives  9$/8r,  9$/90,  92$/9r90. 
Percentage  of  autofrettage. 


INTRODUCTION 


In  this  report  we  present  a simplified  stress  analysis  of  a mandrel 
used  in  the  rotary  forge  process  for  105mm  gun  tubes.  The  objective  was  to 
investigate  the  conditions  of  mandrel  failures  and  the  effects  of 
various  design  parameters  to  enhance  the  life  of  such  mandrels. 

Essentially,  there  are  three  zones  in  the  forging  process,  namely 
sinking,  forging  and  sizing  zones  (Figure  1).  Due  to  physical  con- 
straints, the  sizing  zone  is  subjected  to  higher  stresses1,  and  for 
our  preliminary  analysis,  we  shall  consider  this  zone  under  the  plane 
condition,  i.e.  independent  of  the  axial  dimension. 

To  study  the  crack  propagation,  it  is  necessary  to  consider  the 
most  severe  loading  configuration.  Based  on  the  design  of  the  forging 
process,  it  seems  that  only  two  hammers  are  engaged  with  the  preform 
during  a fraction  of  the  forging  cycle,  causing  peak  stresses  in  the 
mandrel . 

Due  to  this  non-axisymmetric  loading  for  a fraction  of  the  cycle, 
it  is  shown  that  there  are  regions  of  separation  between  the  mandrel 
and  the  preform  causing  high  tensile  stresses  at  the  inner  radius  of 
the  mandrel. 

Assuming  a total  load  transfer  to  the  mandrel  through  the  contact 
region,  the  stresses  are  computed.  Results  compare  favorably  with  a 
finite  element  method. 

1Lahoti,  G.  D.  and  Altan,  R.,  "Analysis  and  Optimization  of  the  Radial 
Forging  Process  for  Manufacturing  Gun  Barrels,"  WVT-CR- 74054,  1974. 
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Figure  1.  A schematic  diagram  of  a mandrel  and  a billet  showing 
sinking,  forging  and  sizing  zones. 
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Thermal  analysis  is  carried  out  based  on  steady  state  conditions 
and  known  heat  flux  of  the  cooling  fluid  through  the  mandrel.  These 
stresses  are  superposed  with  the  mechanical  stresses. 

The  stress  intensity  factors  are  then  computed,  based  on  the  penny 
shaped  crack  under  arbitrary  loading  conditions  and  they  are  compared 
with  two  dimensional  frontal  crack  solutions. 

The  effects  of  change  in  the  internal  diameters  of  the  mandrel 
and  the  effects  of  mandrel  autofrettage  are  studied.  An  optimum 
percentage  of  autofrettage  is  obtained. 

LOSS  OF  CONTACT 

Consider  a smooth  unbonded  inclusion  as  shown  in  Figure  2.  Assume 
that  the  inclusion  has  a higher  elastic  modulus  compared  with  the 
surrounding  unbounded  medium  which  is  subjected  to  compressive  stresses 
as  shown. 


Figure  2.  A mandrel  surrounded  by  an  infinite  medium  under  a compres- 
sive stress,  Sx. 
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Using  the  usual  notations  the  stress  functions  for  the  exterior 
medium  can  be  written  as 

$ = Sxr2  + Sxr2cos20  + Alogr  + Dr  cos20  + Fcos20  (1) 

Using  the  boundary  conditions  of  no  shear  stress  and  no  radial  displace- 
ment at  the  interface  r = b,  we  have 

A = - — (l-2v)b2,  D . ■C1.-2v^'‘Sx  , F = b2Sx 
2 4 (5-6v)  2 (5-6v) 

The  normal  stress  at  the  interface  then  becomes 

ar  = 77  v r {(1-  V)(5-6v)  + 2(4-3v)cos20}  (2) 

[b-oVJ 

From  (2),  it  is  easily  seen  that  normal  stresses  have  a transition 
from  a compression  to  a tension  at  about  0 = 53°  for  v » 1/4.  Since 
the  inclusion  is  unbonded,  this  shows  the  existence  of  separations  at 
the  interface.  To  determine  this  contact  region,  we  have  to  reformu- 
late the  problem  as  done  in  the  sequel. 


2Timoshenko,  S.  and  Goodier,  J.  N.,  THEORY  OF  ELASTICITY,  2nd  Edition, 
McGraw-Hill,  1951. 
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FORMULATION  OF  CONTACT  PROBLEM 


Figure  3.  An  unbonded  inclusion  in  an  infinite  medium  subjected  to 
Sx  showing  regions  of  separation  and  regions  of  contact. 

We  reconsider  the  problem  with  contact  regions  and  separation 

regions  as  shown  in  Figure  3.  The  stress  functions  are  given  by 


$ = - j Sxr2  + j Sxr2cos20  + AQlogr  + 

00  00 

l Anr“ncosn6  + £ Bnr"n+2cosn0  (3) 

n=2,4  n=2,4 

Boundary  conditions  are 

Tr0(b,0)  = 0 O<0<jr/2  (4) 

ur(b,0)  = 0 O<0<n  (5) 

0r(b,0)  = 0 f|<0<7r/2  (6) 
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From  (4) , we  have 


Bn  = " 


(nil)  1 a 
(n-1) 


n>2 


From  boundary  conditions  (6)  and  (5)  using  (7),  we  obtain 


00 

\ Eo  + l Encosn0  = 0,  n>0>7T/2 

2 n=2,4. 


where 


00 

K^Eo  + £ K^Encosn0  = - (l+2cos20),  O<0<ri 

n=2,4. 


a(n+1)An 

bn  , n+2 
b 


1 

4 (1-v) 


k4  . IXS=fi> 1 

2 6 (1-v) 


n+a0 


a„ 


l-2v 

2 (1-v) 


(7) 


(8) 

(9) 


Equations  (8)  and  (9)  are  the  dual  series  relation  for  the  determin- 
ation of  the  unknown  contact  area  as  well  as  unknown  contact  stresses. 
SOLUTION  OF  DUAL  SERIES  AND  THE  CONTACT  ANGLE 

The  dual  series  can  be  solved  by  iteration  procedure  after 
reducing  to  an  integral  equation.  In  this  section  we  outline  the 
solution  procedure  up  to  the  second  iteration.  Using  the  asymptotic 
behavior  of  for  large  n,  the  dual  series  can  be  written  as 


6 


(10) 


1.  EQ  + £ Encosn9  = 0 ti<0£tt/2 

2 n=2,4 


KlEo  + l n Encosn6  = " — (1+2cos20)  l (K^  - i)E  cosnO 

2,4,6  n=2,4  n n (11) 

0<9<n 


Let 


hence 


j E0  + I Encosn0  = H(0)  O<0<n 

2 n=2 ,4 


En  = £ / H(t)cosntdt 


(12) 


(13) 


Substituting  Er  from  (13)  into  (11)  and  by  truncating  the  series  on 
the  right  hand  side  up  to  n = 2,  we  get  the  following  integral  equation 


/ 4 f 1 4 * ‘ cv 

Kx  - / H(t)dt  + - / K(t,0)H(t)dt  ^ — (l+2cos20) 

17  o u o 2 


E2cos20 


Where  the  kernel  is 

K(t,0)  = - log2 | cos20-cos2t | 

The  above  singular  integral  equation  can  be  solved  exactly,  and  the 
procedure  can  be  repeated  for  any  higher  order  approximation.  The 
results  of  the  solution  are  given  in  the  following  table 


(14) 

(15) 


Poisson's  ratio  v 

0.2 

0.3 

0.4 

Contact  angle  q 

54.28° 

49.31° 

44.64° 

The  contact  angle  for  the  mandrel  is  about  50°.  This  is  very  close 
to  the  results  we  obtained  in  the  previous  section  by  a simpler  model. 
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STRESS  DISTRIBUTION  IN  THE  MANDREL 


Because  of  the  dimensions  involved  the  model  analyzed  in  the 
previous  section  can  be  used  for  actual  mandrel-preform  contact.  In 
addition,  it  is  seen  that  the  contact  is  independent  of  the  magnitude 
of  the  load.  Hence,  even  though  more  exact  distribution  of  contact 
stresses  may  be  parabolic,  we  may  consider  a uniform  and  complete  load 
transfer  from  the  hammer  to  the  mandrel,  see  Figure  4.  This  assumption 
can  be  justified  on  the  viscous  flow  model  usually  assumed  in  the  metal 
forming  process'^. 


Figure  4.  A mandrel  is  subjected  to  a uniform  compressive  stress  over 
a part  of  the  outer  boundary. 


3Nadai,  A.,  THEORY  OF  FLOW  AND  FRACTURE  OF  SOLIDS,  2nd  Edition,  McGraw- 
Hill,  1950. 
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With  the  given  dimensions  supplied,  this  stress  h is  computed 
as  38  Ksi  corresponding  to  a load  per  hammer  of  700  tons  with  an 
allowance  of  10%  overload.  Actual  stress  distribution  will  be  non- 
dimensionalized  with  respect  to  this  stress. 

The  stress  function  for  the  mandrel  can  be  written  as 


* ■ Bol°g=>I>0r*  * J ( rn  * A.  r‘"  * -^L  r11*2  * it  r-n*2}cosn0 

n- L*  y T" 


are 

(16) 

Tr0(b,0)  = 0 

O<0  <7T  / 2 

(17) 

Tr0  fa>0^  = 0 

o<0<7 r/2 

(18) 

Or(a,0)  = 0 

0<B<] r/2 

(19) 

-h 

o_(b,0)  = { 
r 0 

o<e<n 

n<0<ir/2 

(20) 

a 1 

From  (17)  and  (18),  with  p = — = — , we  have 

b R 


C„  - -A..-*  - B„a- 

(Rzn-1) 


2n-2  (R2-!) 


R2 (R2n-1) 


(21) 


D = Ana2n-2R2n_2  . . Ra-2  jaf-T-V-b 

(R2n- 1)  R2(R2n-l) 


(22) 


And  from  (19),  we  obtain 


B = [ (n+l)R2n-nR2n~2-llR2a2n 


(R2n+2+(n-l)Rz-n) 


An 


(23) 
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For  convenience,  let  us  define  the  following  quantities 
= - (i-P2n-2)/(i-p2n) 

r2  = p-2(i-p2)/(i-p2n) 

r3  = p2n(p"2-i)/(i-p2n) 

r4  = p-2(i-p2n+2)/(i-p2n) 

r5  = p2n(-p2n-np2+n+l)/[-np2n+2+(n-l)p2n+l] 

r6  = p2^.-^.) 


6 

’7 


1 2 5J 

r7  = p_2n+2(r3-r4r5) 

,-n+2  , 


r8  = -np-n+2(l+r5)-(n-2)p-nr6-(n+2)pnr7 

Using  (21),  (22),  (23)  and  (24),  it  is  seen  that 

00 

^(b,©)  = B ( — »r  _-ir  ) + J ^a  2rscosn0 
° b a n=2,4  8 

From  (20),  using  Fourier  inversion  formula,  we  have 


Bo  = " 


2hri  azb 


2u2 


Ana 


it  (a2-b2) 
n-2  _ 4h  sinnp 


im  r 
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This  formally  completes  the  problem.  The  stress  distribution  of 
interest  to  us  is  the  tangential  stress.  It  is  given  by 

- = 2U  ( 1 + 4)  ♦ l {np_Il+2  (2L)n“2  + 

h it  (1-p  ) p2  r2  n=2,4  n7lT3  b 


(24) 


(25) 


(26) 


np 


-n+2(rrn-2r5  + (n+2)p-n(I.)nr6  + (n-2)pn(^.)-nr7}cosn0 


(27) 


10 


For  the  dimensions  of  the  mandrel,  the  above  stress  distributions  for 
0 = 0 and  0 = tt/2  are  plotted  in  Figure  5.  A similar  finite  element 
model  was  constructed  and  shown  in  Figure  6 and  the  results  for  the 
finite  element  model  are  also  indicated  in  Figure  5. 

THERMAL  STRESSES  IN  THE  MANDREL 

Due  to  lack  of  temperature  data  for  the  mandrel,  several  heat 
transfer  models  have  been  tried  and  ruled  out  because  of  unreasonable 
high  thermal  stresses.  The  final  model  used  in  our  computation  is  a 
simple  one  dimensional,  steady  state  heat  flow  method.  The  temper- 
atures of  the  coolant,  water,  in  and  out  the  mandrel  are  measured,  and 
the  capacity  of  the  water  pump  is  known.  This  gives  the  heat  flux 
taken  away  by  the  coolant  which  is  balanced  by  the  heat  flux  trans- 
ferred from  the  preform  to  the  mandrel.  The  temperature  gradient 
dX 

— at  the  mandrel  preform  interface  is  then  obtained, 
dr 

The  heat  balance  is  given  by 
dT 

k — (2irb)£  = heat  flux  taken  away  by  the  coolant 
Using  measurements  recorded  in  an  actual  manufacturing  run  of  M68  tube 
in  the  rotary  forge:  the  inlet  and  outlet  temperatures  of  cooling 
water  are  65°F  and  80°F,  the  flow  is  80  gallons  per  minute;  taking  the 
thermal  conductivity  k = 0.11  cal/°C-min-sec,  SL  = 368  mm,  we  obtain 

= 3.625  °C/mm  (28) 

r=b 
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r/ b 


Figure  5.  Non-dimensional  stresses  as  a function  of  r/b  for  0=0  and 
tt/2  for  the  mandrel  shown  in  Figure  4. 
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Figure  6.  A finite  element  idealization  of  one  quarter  of  the  mandrel 
shown  in  Figure  4. 
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Assume  a steady  state  condition 

T = c1logr  + c2 

Using  the  boundary  conditions  (28)  and  Ta  = 80°F,  we  have 

Tb-Ta  = 165.27  log(l/p)  , °C 
= 297.49  log(l/p)  , °F 


Defining  parameters  5 and  s as 

/•  \ 


^ ' 


6 = 


Ea(Tb-Ta) 

2 (1-v) log (1/ p) 


s = P2log(l/p) 

.2- 


(1-P  ) 

The  thermal  stress,  normalized  with  respect  to  h,  is  given  by 


0> 


X-  * £ {-l-log(r/b)+s (1  + — 

.*t  h h (r/b)z 

where  6/h  for  the  mandrel  is  1.17. 

TOTAL  STRESS  DISTRIBUTION  AND  STRESS  INTENSITY  FACTORS 

The  total  stress  distribution  for  the  mandrel  is  the  sum  of 
mechanical  and  thermal  stresses.  From  (27)  and  (32),  we  have 

-«$  ♦($ 

n total  1 mechanical  thermal 


(29) 


(30) 


(31) 


(32) 


(33) 


It  should  be  noted  that  this  distribution  of  stress  is  valid  only 
for  a fraction  of  a complete  forging  cycle  when  only  two  hammers  are 
being  engaged.  For  the  rest  of  the  cycle  when  all  four  hammers  are 
engaged  the  stress  distribution  may  not  be  as  severe.  However  there 
exists  a prolonged  thermal  load.  For  the  present  design  the  tangential 
stress  at  9 = 0 (under  the  hammer)  is  plotted  in  Figure  7. 
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Total 


Figure  7.  The  total  tangential  stress  distribution  along  0 = 0 and  its 
linear  approximation  near  the  bore  of  the  mandrel. 


For  the  determination  of  stress  intensity  factors,  we  approximate 
the  stress  distribution  by  a linear  one  in  the  neighborhood  of  the 
inside  radius  as  shown  by  the  dotted  line  in  Figure  7. 

As  far  as  the  stress  intensity  factor  is  concerned,  we  consider  a 
penny  shaped  crack  of  radius  c under  uniform  and  linear  stress  distri- 
butions. For  the  solution  of  the  problem,  we  use  Boussinesq-Papkovitch 
functions  and  reduce  the  problem  to  the  one  function  approach  given 
in  Ref.  [4] . 


where 


2Gur  = (l-2v)$r  + z$rz 
2Gu0  = (1-2V)  I $0  + L $ 

2Guz  = -2 (l-v)$z  + z$zz 
ar  = (l-2v)$rr  + z$zrr  - 2v$ 
a0  = - (l-2v)$rr 


zz 


^zz  + r ^rz  + "pr  ^z00 


a = -$  + zO 

z zz  ^zzz 


(34) 


r0 


ShM 


$ 

r r0 


_ Xl-ftO  $ + - $ _ — $ 

r2  0 r r0z  r2  v0z 


Trz  ” z4rzz 
Tz0  “ T ^0zz 


V2$  = 0 


(35) 


4Green,  A.  E.  and  Zerne,  W. , THEORETICAL  ELASTICITY,  2nd  Edition, 
Oxford,  1968. 


16 


For  the  problem  at  hand,  we  use  Kobayashi  Potentials  as 


oo  . 

$n  = V^110  / K 3 2 Jn+3/2  (Kc)Jn(Kr)e_KZdK 
o 


(36} 


Using  (34)  and  (36)  after  some  manipulation,  it  can  be  shown 


r>c 


2Guz(z=0)  = ^ 


(c2-r2)1/2  rn 
2 (1"V;i  r-  n+3/2  3 kncosnB  > rlc 

/2  c •5/zr(-) 

2 


(37) 


pr  n-1/2 , 2 2-v-1/2. 

✓2  c (r  -c11)  _ An 


a7(z=0)  =< 


Up  - 1 ^ 

r r(-) 


f2n+llcn 

^2?(n*  I)rn+1 


F(ji  n+i,  , 

r>c 

(38) 


(2n+l)T (n+1/2) 

n+3/I 


Tl 


n 

r cosn6 


r<c 


T (n+l)/2  c 

The  reason  for  selecting  the  potential  given  in  (36)  becomes  clear  if 
we  study  (37)  and  (38) . It  is  seen  that  the  stress  has  the  crack 
singularity  at  r = c and  the  displacement  vanishes  outside  of  r = c. 
These  are  precisely  the  conditions  required  for  a penny  shaped  crack 
under  normal  stress  of  the  form  rncosn0. 

For  our  purpose,  it  is  sufficient  to  consider  only  the  two-term 
approximation,  i.e.  $ = where  $q,  are  given  by  n = 0,1  in  (36); 

az(z=0)  = a+e£=a  + 3(|-)cos0  (39) 


^Sneddon,  I.  N.,  MIXED  BOUNDARY  VALUE  PROBLEMS  IN  POTENTIAL  THEORY, 
John  Wiley  § Sons,  1966. 
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Using  (38),  we  have 


From  the  definition  of  stress  intensity  factors,  we  have 


(40) 


Kx  = — (Aq+Aj)  (41) 

It  is  known  that  for  a semi-circular  surface  flaw,  the  stress  intensity 
factor  at  the  deepest  end  of  the  crack  very  closely  approximates  the 
stress  intensity  factor  for  penny  shaped  crack  in  an  unbounded  medium 
provided  the  load  distribution  is  similar.  Hence,  to  a very  good 
approximation,  we  take  Oz  in  equation  (39)  to  represent  a linear 
approximation  of  (ag)  ^ (33).  The  result  for  the  present  mandrel 

design  (a  = 0.9845”)  of  such  computations  is  shown  in  Figure  8.  It  is 
seen  that  the  critical  crack  size  of  the  semi-circular  shape  is  about 
20  thousandths  of  an  inch  for  a fracture  toughness  of  30  Ksi/In.  Even 
though  the  cracks  at  failure  seem  to  be  of  the  semi-circular  shape, 
we  computed  the  stress  intensity  factor  for  a frontal  crack  by  a 
similar  analysis,  i.e. 

K1  ■ ^ <°iFi<P7>  * Wb??1  <42> 

where  and  a 2 correspond  to  uniform  tension  and  bending  stresses, 
respectively,  in  the  linear  approximation  of  (0q) total , F1  F2 
are  given  by  Ref.  [6], 


^Tada,  H.,  Paris,  P.,  and  Irwin,  G. , THE  STRESS  ANALYSIS  OF  CRACKS 
HANDBOOK,  Del  Research  Corp.,  1973. 


18 


(wf/sx) 


Figure  8.  Stress  intensity  factors  as  a function  of  c for  various  values 
of  the  bore  radius  a of  the  mandrel. 
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Fjfx)  = 1.12  - 0.231x  + 10.55x2  - 21.72x3  + 30.39x4 


(43) 


F2(x)  = 1.122  - 1.4x  + 7. 33x2  - 13.08x3  + 14. Ox4 
The  results  are  shown  in  Figure  9.  It  can  be  seen  from  Figures  8 and 
9 that  the  stress  intensity  factors  are  much  higher  for  frontal  cracks 
for  the  same  c. 

EFFECT  OF  VARIATION  OF  THE  INNER  RADIUS  OF  THE  MANDREL 

Similar  analysis  was  carried  out  by  fixing  the  outer  radius  of 
the  mandrel  and  varying  the  inner  radius.  The  outer  radius  was  fixed 
due  to  the  design  constraint,  and  the  total  heat  flux  was  assumed  the 
same  even  with  a slight  reduction  in  the  inner  radius. 

The  total  stress  at  0 = 0 was  computed  for  each  value  of  inner 
radii  and  the  non-dimensionalized  results  are  shown  in  Figure  10.  It 
is  seen  that  the  maximum  tensile  stress  at  the  bore  can  be  reduced  by 
30%  for  the  inner  diameter  a = 1.4  inch.  The  corresponding  stress 
intensity  factors  are  shown  in  figure  8.  It  can  be  seen  that  the 
critical  crack  size  is  now  increased  by  150%  for  the  same  inner  radius 
of  a * 1.4". 

At  this  point,  it  is  worhtwhile  to  compare  the  number  of  remaining 
cycles  the  mandrel  can  sustain  this  load  without  a catastrophic  fracture. 
Assuming  an  undetactable  crack  of  0.003"  and  using  the  propagation 
equation  for  gun  steel,  approximate  results  are  given  below.  In  the 
table,  1600  cycles  is  used  to  obtain  the  number  of  tubes  from  the 
number  of  cycles. 
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STRESS  JKTEUS/TY  FACTOR  R?  ( K si J7* ) 


CRACK  bcprR  C (in.) 


Figure  9.  Stress  intensity  factors  as  a function  of  the  depth 
c of  a frontal  crack. 
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Figure  10.  The  total  tangential  stress  distribution  along  0=0*  for 
various  bore  sizes  of  the  mandrel. 


Inner  Radius 
a inch 

Remaining 

Life  to  Failure 

Ratio  of 
Mandrel  Life 

In  Cycles 

In  No.  of  Tubes 

0.9845 

3530 

2 

1 

0.9 

6580 

4 

2 

0.8 

10880 

7 

3.5 

0.7 

13350 

8 

4 

This  table  shows  that  by  changing  the  inner  diameter  from  1.969" 
to  1.4",  there  is  a four- fold  increase  in  the  life  of  the  mandrel. 
Actual  life  may  be  longer  than  those  indicated  in  the  table  because 
of  our  assumption  of  worst  conditions  during  each  cycle. 

EFFECT  OF  MANDREL  AUTOFRETTAGE 


An  examination  of  Figure  7 indicates  that  the  stress  distribution 
is  exactly  the  opposite  of  the  residual  stress  distribution  of  a fully 
autofrettaged  tube.  This  suggests  the  beneficial  effect  of  having 
autofrettaged  mandrels. 

Let  a)  be  the  percentage  of  autofrettage,  the  residual  stress  is 


O 

y t co 

'Vresidue  = JS  {21°g  Z * 1 * b7  ' 


2 2 
£1  a2 

b2  r? 


1 - 


(J  OJ 

2 (2iog  7 + i - rr)> 


b^ 


2 

a bF' 

a<r<(j  (44) 


Thus  from  (44)  and  (33),  we  obtain  the  total  stresses  and  the  results 
are  plotted  in  Figure  11,  showing  a dramatic  reduction  in  tensile 
stresses. 

The  corresponding  stress  intensity  factors  are  shown  in  Figure  12, 
again  showing  the  dramatic  decrease.  As  done  previously,  we  compute 
the  remaining  life  to  failure  and  the  results  are  tabulated  in  the 
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Figure  11.  The  total  tangential  stress  distribution  along  9=0  for  the 

mandrel  of  present  design  but  with  various  percentages  of 
autofrettage. 
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Figure  12. 
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Stress  intensity  factors  as  a function  of  c for  the  mandrel  of 
present  design  but  with  various  percentages  of  autofrettage. 


40%  autofrettage 


following  table. 


Percentage  of 
Autofrettage  oj 

Remaining 

Life  to  Failure 

Ratio  of 
Mandrel  Life 

In  Cycles 

In  No.  of  Tubes 

0% 

3530 

2 

i 

40% 

47770 

29 

15 

50% 

154000 

96 

48 

60% 

218000 

136 

68 

80% 

885000 

554 

277 

100% 

1660000 

1040 

520 

The  preceding  table  indicates  a dramatic  increase  in  the  mandrel 
lifs.  However,  this  was  based  on  the  analysis  valid  only  for  a fraction 
of  a forging  cycle.  The  rest  of  the  cycle  induces  a large  compressive 
stress  at  the  inner  bore.  The  sum  of  this  compressive  stress  and  the 
compressive  stress  due  to  atuofrettage  may  cause  a reverse  yielding 
phenomenon  at  the  inner  bore  of  the  mandrel.  Such  an  analysis  was 
carried  out  and  it  was  found  that  50%  autofrettage  gives  us  the 
optimum  value. 

In  Figure  13,  we  also  plotted  the  stress  distribution  for  a 
solid  mandrel.  This  shows  negligible  stresses  in  the  mandrel. 

DISCUSSION 

In  this  report  we  have  done  only  a preliminary  analysis  in  two 
dimensions.  The  proper  analysis  should  include  some  three-dimensional 
effects  and  some  non-linear  effects.  However,  before  such  an  extensive 
analysis  is  undertaken,  it  is  necessary  to  apply  proper  boundary  and 
working  conditions.  For  example,  proper  temperature  distribution 
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Figure  13.  The  tangential  stresses  as  a function  of  r/b  for  0=0  and 
■jt/2  for  a solid  mandrel  subjected  to  a loading  shown  in 
Figure  4 with  n = 50°. 
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should  be  measured,  exact  amount  of  heat  removed  should  be  determined, 
and  proper  calibration  of  load  per  hammer,  etc.,  should  be  given. 
CONCLUSION 

It  is  seen  that  by  changing  the  inside  diameter  to  1.4",  we 
increase  the  life  of  the  mandrel  by  a factor  of  four.  On  the  other 
hand,  40-50%  autofrettage  increases  the  life  of  the  mandrel  by  a 
factor  of  15-45.  The  computations  were  based  on  an  initial  crack 
length  of  0.003".  To  further  increase  the  life,  cleaner  material  with 
higher  fracture  toughness  should  be  used. 
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